Molecular Elasticity and the Geometric Phase 
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We present a method for solving the Worm Like Chain (WLC) model for twisting semiflexible 
polymers to any desired accuracy. We show that the WLC free energy is a periodic function of the 
applied twist with period 4-k. We develop an analogy between WLC elasticity and the geometric 
phase of a spin i system. These analogies are used to predict elastic properties of twist-storing 
polymers. We graphically display the elastic response of a single molecule to an applied torque. 
This study is relevant to mechanical properties of biopolymers like DNA. 



PACS numbers: 82.37.-j, 03.65. Vf,87.14.-G,87.15.-V 

Some long molecules are as stiff as needles, others 
as flexible as thread. The elasticity of a stiff molecule 
is dominated by energy, while the elasticity of a flexi- 
ble molecule is dominated by its conflgurational entropy. 
It has lately become feasible [1] to stretch, bend and 
twist single molecules to study their elastic properties. 
The subject of this paper is the elasticity of semiflexible 
polymers in which there is competition between ener- 
getic and entropic effects. We consider a polymer which 
can bend as well as twist [2]. The flexibility of such a 
molecule is characterized by two dimensionless param- 
eters a — Lbp/Ltp and [3 — L/Lbp, where L is the 
length of the molecule and Lbp and Ltp are the bend 
and twist persistence lengths. For example DNA has a 
bend persistence length of about 53nm and a twist per- 
sistence length of about 70nm. 

The main purpose of this letter is to draw attention 
to an experimentally relevant topological subtlety which 
has not been discussed in previous theoretical treatments 
of the elasticity of twisting polymers. To appreciate the 
point, take a strip of paper (a ribbon or belt will do as 
well) and tape one end of the strip to a table. Pull the 
strip taut by its other end and twist it by four half-turns 
(47r rotation) . If you now slacken the strip, you will find 
that it is possible, keeping the end fixed, to pass the 
strip around the end. Pulling the strip taut will reveal 
that the An twist has been released. This demonstration 
shows that a polymer can release twist, two turns at a 
time by going around its end. A twist of 27r cannot be so 
released but can be transformed to — 27r. This of course, 
merely illustrates the well known mathematical fact that 
the rotation group is doubly connected [3]. We will now 
see that this mathematical fact has concrete experimen- 
tal consequences for the elastic properties of twist-storing 
polymers. Mentally replace the table with a translation 
stage, the strip with DNA molecule with a micron sized 
magnetic bead at its free end and let the twisting be done 
with a magnetic field. If the DNA molecule is about 16/i 
long, and not pulled taut, it can release 47r worth of twist 
by passing around the bead. One will find that twisting 
the bead by 47r is equivalent to not twisting it at all! Un- 
der the influence of thermal agitation, the molecule will 



explore all conflgurations which it can reach by contin- 
uous deformation. As we have shown, a configuration 
with a 47r twist can be continuously deformed to a con- 
figuration with twist. It follows that the free energies 
for these two situations are the same. If one were to mea- 
sure the free energy G{il') of the molecule as a function 
of applied twist ip (by say measuring the torque-twist re- 
lation), one would find that G{Tp) is a periodic function 
of ip with period 47r. From the 47r periodicity of the free 
energy, it follows that the torque-twist relation and other 
measurable elastic properties have the same periodicity. 

The above discussion is even more relevant to theo- 
retical models of polymers which do not incorporate self 
avoidance. Such a polymer is a 'phantom chain' [4] and 
can pass through itself. In the absence of self avoidance, 
a polymer does not even need to pass around its end: it 
can pass through itself and so release twist two turns at a 
time. Indeed, such effects have been seen in recent single 
molecule experiments [5] and even earlier [6] : the enzyme 
topoisomerase II converts a real DNA into a 'phantom 
chain' [5,7] and in the presence of this enzyme (which 
plays a crucial role in replication), the DNA molecule re- 
leases twist two turns at a time [5,6]. The release of twist 
through bending modes (geometric untwisting) has also 
been discussed in [8] . 

It is thus clear that in theoretical models which do not 
have self avoidance, the free energy is a periodic func- 
tion of applied twist with period 47r. Bearing this in 
mind, we will now explore the Worm Like Chain (WLC) 
[9] model for twisting polymers, which ignores self avoid- 
ance. As was emphasized recently [10] in a review of 
single molecule experiments, "the precision of control 
and quantitative measurement and simple interpretation 
of these experiments make detailed theoretical analyses 
appropriate" . While there has been some theoretical 
progress [4,11] on the WLC for twisting polymers, in- 
terest has been confined primarily to the high tension 
regime, which is theoretically more tractable because the 
molecule has only small perturbations about a linear con- 
figuration. Experiments of Bensimon [12] have also ex- 
plored this regime. In contrast the present letter theo- 
retically explores the non linear regime of small forces. 
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where the pcrturbativc methods described in [4,11] are 
inappUcable. In this regime we find that the WLC free 
energy has a iir periodicity in contrast to the aperiodic 
free energy of the high tension regime. 

We first describe the WLC model and summarize the 
existing theoretical analyses of this model. We then point 
out a physically relevant mathematical subtlety without 
which the solution of the WLC model is incorrect. We 
then correctly solve the WLC model by a combination of 
analytical and numerical techniques and elucidate some 
of the elastic properties that emerge from the model in 
graphical form. A concluding discussion interprets the 
results. 

WLC Model: The WLC model ignores self avoid- 
ance and views the polymer as a framed space curve 
C ~ {x{s) , ei{s)} , i = 1, 2, 3, < s < L of contour length 
L, with an energy cost £(C) for bending and twisting. 
We suppose that one end of the polymer is tethered to 
the origin x(0) = and the other end at x{L) = r is 
tagged. The unit tangent vector 63 = dx/ds to the curve 
describes the bending of the polymer while the twisting 
is captured by a unit vector ei normal to 63. 62 is then 
fixed by 62 = 63 x ei to complete the right handed mov- 
ing frame ei{s),i = 1,2,3. The rate of change of the 
moving frame ei{s) along the curve can be measured by 
its "angular velocity vector" Q defined by 

-^ei(s) = n X ei{s). (1) 
as 

The components of in the moving frame are fi^ = fl.Ci 
and the energy £{C) of a configuration C is given by 

£[C] = 1/2 [\A{{nif + {Q^f) + C{nsf]ds, (2) 
Jo 

where A is the bending modulus and C the twist modu- 
lus. 

Imagine that the ends of the polymer and the frames 
at these ends (ei(0), ei(L)) are held fixed. We wish to 
compute the number of configurations 

Q(r, 6^(0), ei{L)) = Sc exp(-f [C]/fcT), (3) 

counted with Boltzmann weight cxp{—£[C]/kT), which 
start at the origin with initial frame 6^(0) and end at 
r with final frame ei{L). The function Q is related to 
the free energy of the molecule and its measurable elas- 
tic properties like the force-extension relation (FER) and 
the torque- twist relation (TTR). An overall multiplica- 
tive constant is not important in the calculation of Q. 
This only leads to an additive constant in the free en- 
ergy, which drops out on differentiation and does not 
affect elastic properties. 

Let us fix a "lab frame" and write ef{s) = R'^b{s)ei, 
where a = 1,2,3 is a vector index (as opposed to the 



frame index i) and R{s) E SO{3) is a 3 x 3 rotation ma- 
trix. There is a clear analogy between the elastic prop- 
erties of the WLC and the motion of a top. Indeed as 
Mezard and Bouchiat [11] and Moroz and Nelson [4] point 
out, the WLC problem can be mapped to the quantum 
mechanics of a symmetric top. From this mapping, one 
may naively conclude [13] that the periodicity of the free 
energy is 2n. We now show that a careful treatment of 
the path integral (3) gives the correct 47r periodicity. 

Each configuration C in (3) is characterized by a curve 
{i?(s)} in the rotation group with fixed end points -R(O) 
and R{L). The sum is over all configurations which are 
sampled under the influence of thermal agitation. It is ev- 
ident that we should only sum over configurations {i?(s)} 
in a single homotopy class: thermal agitation only causes 
continuous deformations and therefore cannot knock the 
polymer out of its homotopy class. Q thus depends not 
only on (f, -R(O), R{L)) , but also on the homotopy class 
[{i?(s)}] of the paths {R{s)} being summed over. This 
information is exactly captured by going to the cover- 
ing space SU{2) of the rotation group SO{3). This step 
is essential to correctly describe the elastic properties of 
the twisting polymer. (If we do not take this step but 
remain on S0{3), we are effectively summing over both 
homotopy classes, which is a physically incorrect proce- 
dure.) The result is that while the WLC Hamiltonian 
is the same as that of the top, the WLC configuration 
space is not the configuration space 5*0(3) of the top, 
but its double cover SU{2) [3]. As we will see below, 
this results in a Air periodicity for the free energy. SU (2) 
is the same as S*^, the four dimensional sphere defined 
by = l,2,3,4,},Sa(a;")2 = 1. In fact, are 

the Cayley-Klein parameters [14] traditionally used in 
describing tops. 

Let {g(s)} be a continuous curve in SU{2) which maps 
down to the curve {R{s)}. g{s) satisfies the differential 
equation '^^^ = i/2{Q..<j)g{s)^ whose solution is a path 
ordered exponential: g{s) = P[exp iQ{s').a /2ds']g{0). 
C is now described by a curve {g{s} in SU{2), with fixed 
end points g(0) and g{L). The standard Euler angles 
{9, (j), tp) on the rotation group can be used as co-ordinates 
on SU{2) = if the range of tp is extended to Att [3]. 
SU{2) acts on itself by right and left action generated 
[15] by space fixed {J^, Jy, Jz) and body fixed (Ji, J2, J3) 
angular momenta. 

We can now write (3) more correctly as Q(?^, <Zo, 'Zl) 
where go = 5(0) and ql = g{L), to explicitly display the 
homotopy class dependence of Q. Q{r,qo,qL) has the 
path integral representation: 

J D[<7(s)]e[-^(^)/'=^^l(5(f(L) - f). (4) 

TV is a normalisation constant and the path integral is 
over all paths that go from qo to qL on S^. We now 
pass from Q{f,qQ,qL) to its Laplace transform defined 
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as Q{f,qo,qL) = J d'rexp[f.r/LBp]Q{r,qo,qL), where 
Lbp = A/kT. Performing the elementary integrations 
and changing variables to r = s/ Lbp and co = QLbp, we 
sec that Q{f, qa, qh) can be represented as MZ[f, Qq, Ql), 
where Z has the path integral representation 

where a = A/C = Lbp/Ltp- This is clearly the quan- 
tum amplitude < qL\exp[—PHf]\q{) > for a particle on 
the surface of a 3-sphere to go from an initial position go 
on to a final position gi, in imaginary time /3 in the 
presence of an external force field. The Hamiltonian is 
Hf ^ Ho- fcose, where Hq = l/2{Ji^ + J^^ + aJ^), 
which is the Hamiltonian of a symmetric top. If the exact 
eigenstates of B.^ were known, we could write: 

Z = S„ exp[-/JE„]u* (go)wn(gL), (6) 

where {^/^(g')} is a complete set of normalised eigenstates 
of the Hamiltonian Hj and are the corresponding 
eigenvalues. Even though Hj cannot be diagonalised an- 
alytically, we can exploit its symmetries to reduce the 
problem to a numerically tractable form. and J3 com- 
mute with the Hamiltonian iJj, reflecting the symme- 
try under space-fixed and body fixed rotations about the 
third axis. As a result, Z depends only on the differ- 
ences (/) = (f>L — 4'o and tp = tpL — fpo and we write 
Z{f,9o,9L,(f>,'ip). Consider the dependence of Z on tp. 
Since all the wave functions in (6) are single valued func- 
tions on S^, it follows that Z{tp) is periodic in ip with 
period Aw: Z{ip + 4tt) = Z{tjj). This means that the free 
energy G = —1/ (3Log[Z] and all elastic properties have 
the same period. This is the first main result of this let- 
ter. The fact that the periodicity is Air and not 27r can be 
traced to the fact the sum in (6) extends not only over 
tensorial states but also over spinorial ones. The vari- 
ation of G with respect to the variables [f ,9o, 9^,(1), 
gives the elastic response to stretch (/), bend {9^(1)) and 
twist (V'). G can be computed numerically for any value 
of its arguments using mathematica programs [16] that 
run for a few minutes on a PC, using methods similar to 
those of [17]. From G we can extract all possible infor- 
mation regarding the elasticity of a polymer with bend 
and twist degrees of freedom. For instance one can pre- 
dict the form of extension versus twist curves for various 
values of the stretching force. 

The 47r periodicity of the free energy strongly moti- 
vates the use of spinorial methods. In fact, the WLC 
configuration space is the same as the set of nor- 
malised states of a spin- 1/2 quantum system. We will 
now show that there is a mapping between a configu- 
ration of a twisting polymer and the quantum evolu- 
tion of a spin- 1/2 system. This brings out an inter- 
esting connection between WLC elasticity and the Ge- 
ometric Phase. Let us introduce a 2 component com- 
plex vector (a spinor) = + ix^,^^ = + ix"^, 



which is normalised (^^^ = 1). We can write = 
cos9/2exp-"f'/^ exp*^/^ = sin9/2exp"f'/^ exp*'^/^ a,nd 
thus introduce co-ordinates {0,<p,'4i) ranging from to 
(7r,27r,47r) respectively. These are similar to Euler an- 
gles on the rotation group and differ only in the range 
of ip. The frame can be expressed as 63 = ^^<t^, 
ei -|- ie2 = ^"^(io-2)(7^, where the crs arc the usual Pauli 
matrices. Notice that altering tp by 27r flips only the 
sign of ^ and therefore does not affect the frame. Using 
this mapping between 2 component spinors and frames, 
we can import ideas from the geometric phase to under- 
stand WLC elasticity. The information in the spinor ^ 
can be decomposed into an overall phase ^p/2 describing 
twist and a ray 63 describing bend. Fix a conflguration 
C and note that ^(s) = g(s)^(0) satisfles the Schrodinger 
differential equation i^^^j^ = ^^(s), where h = —i}.a/2 
is the "Hamiltonian" of a spin half particle in an external 
magnetic fleld Q. We can now decompose the difference 
■0/2 = ip{L)/2 — ■0(O)/2 between the final and initial 
phases into a geometric phase and a dynamical phase. 
The dynamical phase is given by the integral of the ex- 
pectation value of the "Hamiltonian" dsS,'^h^, which, 

using the definition of 63 above is seen to be dsfl3/2, 
half the twist il^ integrated along the polymer. The ge- 
ometric phase is given by half the solid angle swept out 
by the ray e^{s). If the initial and final rays are distinct 
(but not antipodal), one can join them by the unique 
shorter geodesic [18] to enclose a solid angle. The total 
twist difference ijj is the sum of the "dynamical twist" -the 
integrated local twist- J fl^ds and the "geometric twist"- 
the solid angle swept out by the tangent vector- which 
depends on the bending of the polymer. In the literature 
[4,11,19], the distribution of applied twist between twist- 
ing and bending is compared with the decomposition of 
link into twist and writhe. This result is referred to as 
White's theorem [20] (though an earlier reference is [21]). 
The discussion [4,11,19] applies to closed, self-avoiding 
polymers which have been twisted an integral number of 
times. In contrast, our treatment applies also to open 
polymers which have been twisted a fractional number 
of times. However, since the WLC model does not take 
into account self avoidance, a twist of An is equivalent 
to no twist and the integral part of the twist is only 
measured modulo two. Our treatment captures the frac- 
tional part of the applied twist (which is geometrical) and 
the earlier treatment captures its integral part (which is 
topological). In this sense, the two discussions are com- 
plementary. The analogy between twist elasticity and 
the geometric phase is the second main result of this let- 
ter. The analogy has also been noted in Ref [19] , which, 
however, uses a vectorial correspondence rather than a 
spinorial one. The decomposition of applied twist into 
a geometrical and a dynamical part leads to a coupling 
between the bend and the twist degrees of freedom and 
has a direct bearing on the elastic properties of the WLC 
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polymer. As a specific illustration we give the results for 
the special case of pure twist elasticity. 

Pure twist elasticity: We suppose that the tagged end 
is not constrained in position, but only in orientation. 
Integrating Q{r,qo,<lL) over r, we see that the applied 
force / vanishes. We also suppose that the initial 63(0) 
and final es(L) tangent vectors are both in the same di- 
rection (which we take to be the z direction). We com- 
pute the distribution Z{ip) of ip. In this case only states 
for which m = g contribute [22] and Z takes the form 
Z = YigZgg. Uslug staudard techniques from angular 
momentum theory, we can express Z bs Z = TigC^^^Qg 
where Qg = Ef^^g0 + l)exp[-f3/2ij{j + l) + {a-l)g% 
where j runs in integer steps. By an inverse Fourier trans- 
form, we compute P{tjj) and the free energy G{ip) = 
— l/f3Log[P{tp)]. By differentiating with respect to tp, 
we compute the torque Q = dF/dtp needed to twist the 
molecule by an angle ip. The torque-twist relation is plot- 
ted in Fig.l. This graph, which is the third m,ain result 
of this letter, describes the pure twist elastic properties 
of a molecule in the WLC model. These graphs are eas- 
ily interpreted in terms of the geometric phase ideas de- 
scribed earlier. For large a {a is the ratio of the bend 
to the twist persistence length), twist costs very little 
energy, the molecule twists without bending, and as it 
takes hardly any torque to twist the molecule, the TTR 
is almost flat. As a decreases, the applied twist is shared 
between the twist and the bend. When a is zero twist- 
ing is prohibitively expensive and the applied twist is all 
taken up by the bend. This causes the molecule to buckle 
just as a towel does when it is wrung. When a = 0, 
the "polarisation vector" ei is parallel transported along 
the polymer. The distribution Z{'4!) then reduces to the 
distribution of solid angles (Berry Phases) enclosed by 
closed Brownian paths on the (Poincare) sphere, which 
was calculated in [23] in the context of depolarised light 
scattering. 

In this letter we have solved the WLC model with 
bend and twist degrees of freedom and noticed analo- 
gies to spinl/2 systems and the geometric phase. These 
analogies lead to a description in terms of a particle on a 
sphere in external gravitational and magnetic fields. Such 
analogies, apart from giving us analytic tools to solve the 
problem virtually exactly for the first time, also provide 
simple physical pictures: Imposing a twist on a molecule 
is like applying a magnetic field. The helical shape of a 
towel when it is wrung is similar to the helical trajectory 
of a particle in a magnetic field. We hope this letter will 
encourage experimental work on twisting polymers in the 
nonlinear low tension regime and set up a dialog between 
the theory and experiments on molecular elasticity. 
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FIG. 1. The Torque-Twist relation for (3 = L/L 
= L/Ltp = 0, 1, 7. 



